Phase gate of one qubit simultaneously controlling n qubits in a cavity 
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We propose how to realize a three-step controUed-phase gate of one qubit simultaneously control- 
ling n qubits in a cavity or coupled to a resonator. The n two-qubit controUed-phase gates, forming 
this multiqubit phase gate, can be performed simultaneously. The operation time of this phase gate 
is independent of the number n of qubits. This phase gate controlling at once n qubits is insensi- 
^ . five to the initial state of the cavity mode and can be used to produce an analogous CNOT gate 

O ' simultaneously acting on n qubits. We present two alternative approaches to implement this gate. 

' One approach is based on tuning the qubit frequency while the other method tunes the resonator 

frequency. Using superconducting qubits coupled to a resonator as an example, we show how to 
00 ' implement the proposed gate with one superconducting qubit simultaneously controlling n qubits 

selected from A'^ qubits coupled to a resonator (1 < n < A^). We also give a discussion on realizing 
' (— I ', the proposed gate with atoms, by using one cavity initially in an arbitrary state. 
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I. INTRODUCTION 



Quantum information processing has attracted considerable interest during the past decade. The building blocks 
of quantum computing are single-qubit and two-qubit logic gates. So far, a large number of theoretical proposals 
' for reahzing two-qubit gates in many physical systems, such as trapped ions, atoms in cavity QED, nuclear mag- 
, netic resonance (NMR), and solid-state devices, have been proposed. Moreover, two-qubit controlled- not (CNOT), 
controUed-phase (CP), zSWAP gates, or other two-qubit operations have been experimentally demonstrated in ion 
, traps [1], NMR [2], quantum dots [3], and superconducting qubits [4-8]. 

• In recent years, analogs of cavity quantum electrodynamics (QED) have been studied using superconducting Joseph- 
son devices. There, a superconducting resonator provides a quantized cavity field, superconducting qubits act as 
I artificial atoms, and the strong coupling between the field and superconducting qubits has been demonstrated [9]. 
' Theoretically, many approaches for performing two-qubit logic gates using superconducting qubits coupled to a super- 
: conducting resonator or cavity have been proposed [e.g., 10-15]. Moveover, experimental demonstrations of two-qubit 
gates or two-qubit operations using superconducting qubits coupled to cavities have recently been reported [7,8]. 
. Attention is now shifting to the physical realization of multi-qubit gates (see, e.g., Ref. [16]) instead of just two-qabii 

^ ■ gates. It is known that any multi-qubit gate can be decomposed into two-qubit gates and one-qubit gates. When 
J-j using the conventional gate-decomposition protocols to construct a multi-qubit gate [17,18], the procedure usually 
. . . I becomes complicated as the number of qubits increases. Therefore, building a multi-qubit gate may become very 
difficult since each elementary gate requires turning on and off a given Hamiltonian for a certain period of time, and 
each additional basic gate adds experimental complications and the possibility of more errors. During the past few 
years, several methods for constructing multi-qubit phase gates with n-control qubits acting on one target qubit based 
on ion traps [19], cavity QED [20,21], or circuit QED [22] have been proposed. As is well known, a multi-qubit gate 
with multiple control qubits controlling a single qubit plays a significant role in quantum information processing, such 
as quantum algorithms [e.g. ,23,24] and quantum error-correction protocols [25]. 

In this work, we focus on another type of multi-qubit gates, i.e., a multi-qubit phase gate with one control qubit 
simultaneously controlling n target qubits. This multi-qubit gate is useful in quantum information processing such as 
entanglement preparation [26], error correction [27], quantum algorithms (e.g., the Discrete Cosine Transform [28]), 
and quantum cloning [29]. In the following, we will propose a way for realizing this multi-qubit gate using (n + 1) 
qubits in a cavity or coupled to a resonator. To implement this gate, we construct an effective Hamiltonian which 
contains interaction terms between the control qubit and each subordinate or target qubit. We will denote this n- 
target-qubit control-phase gate as an NTCP gate [Fig. 1(a)]. We present two alternative approaches to implement 
this NTCP gate. One approach is based on tuning the qubit frequency, while the other method tunes the cavity (or 
resonator) frequency. We present these two alternative methods because some experimental implementations might 
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FIG. 1: (Color online) (a) Schematic circuit of a n-target-qubit control phase (NTCP) gate with qubit 1 simultaneously 
controlling n target qubits (2, 3, n+1). The NTCP gate is equivalent to n two-qubit control phase (CP) gates each having a 
shared control qubit (qubit 1) but a different target qubit (either qubit 2, 3, or n+ 1). Here, Z represents a controUed-phase 
flip on each target qubit, with respect to the computational basis formed by the two eigenstates |+) and |— ) of the Pauli 
operator ax- Namely, if the control qubit 1 is in the state |— ), then the state |— ) at each Z is phase-flipped as |— ) — )■ — ]— ), 
while the state |+) remains unchanged, (b) Relationship between a n-target-qubit controUed-NOT gate and a NTCP gate. 
The circuit on the left side of (b) is equivalent to the circuit on the right side of (b). For the circuit on the left side, the symbol 
© represents a NOT gate on each target qubit, with respect to the computational basis formed by the two eigenstates |0) and 
|1) of the Pauli operator a^- If the control qubit 1 is in the state then the state at © is bit flipped as |1) — ^ |0) and |0) 
— > However, when the control qubit 1 is in the state |0), the state at © remains unchanged. On the other hand, for the 
circuit on the right side, the part enclosed in the (red) dashed-line box represents a NTCP gate. The element containing H 
corresponds to a Hadamard transformation described by |0) — >■ |+) — (l/\/2) (|0) + |1)), and |1) — > |— ) — (1/a/2) (|0) — |1)) . 



find it easier to tune the qubit frequency, while others might prefer to tune the cavity frequency. For solid state 
qubits such as superconducting qubits and semiconductor quantum dots, the qubit transition frequency (or the qubit 
level spacings) can be readily adjusted by varying the external parameters [30-34] (e.g., the external magnetic flux 
for superconducting charge qubits, the flux bias or current bias in the case of superconducting phase qubits and flux 
qubits, see e.g. [30-33], and the external electric field for semiconductor quantum dots [34]). Also, the cavity mode 
frequency can be changed in various experiments (e.g., [35-39]). 

As shown below, our proposal has the following advantages: (i) The n two-qubit CP gates involved in the NTCP 
gate can be performed simultaneously; (ii) The operation time required for the gate implementation is independent of 
the number n of qubits; (iii) This proposal is insensitive to the initial state of the cavity mode, and thus no preparation 
for the initial state of the cavity mode is needed; (iv) No measurement on the qubits or the cavity mode is needed 
and thus the operation is simplified; and (v) The proposal requires only three steps of operations. 

Note that a CNOT gate of one qubit simultaneously controlling n qubits, shown in Fig. 1(b), can also be achieved 
using the present proposal. This is because the n-target-qubit CNOT gate is equivalent to an NTCP gate plus two 
Hadamard gates on the control qubit [Fig. 1(b)]. 

To the best of our knowledge, our proposal is the only one so far to demonstrate that a powerful phase gate, 
synchronously controlling n qubits, can be achieved in a cavity or resonator, which can be initially in an arbitrary 
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state. This proposal is quite general and can be applied to physical systems such as trapped atoms, quantum dots, 
and superconducting qubits. We believe that this work is of general interest and significance because it provides a 
protocol for performing a controlled-phase (or controUed-not) gate with multiple target qubits. 

This paper is organized as follows. In Sec. II, we introduce the n-target-qubit control phase gate studied in this 
work. In Sec. Ill, we discuss how to obtain the time-evolution operators for a qubit system interacting with a single 
cavity mode and driven by a classical pulse. In Sec. IV, using the time-evolution operators obtained, we present two 
alternative approaches for realizing the NTCP gate with qubits in a cavity or coupled to a resonator. In the appendix, 
we provide guidelines on how to protect multi-level qubits from leaking out of the computational subspace. In Sec. V, 
using superconducting qubits coupled to a resonator, we show how to apply our general proposal to implement the 
proposed NTCP gate, and then discuss its feasibility based on current experiments in superconducting quantum 
circuits. In Sec. VI, we discuss how to extend the present proposal to implement the NTCP gate with trapped atomic 
qubits, by using one cavity. A concluding summary is provided in Sec. VII. 

II. ONE QUBIT SIMULTANEOUSLY CONTROLLING 71 TARGET QUBITS 

For two qubits, there are a total of four computational basis states, denoted by |-|— 1-) , |H — ) , | — h) , and | — ) . 
Here, |-|-) and |— ) are the two eigenstates of the Pauli operator ax- A two-qubit CP gate is defined as follows 

^ |+-)^|+-), 

^ I— )^-|— ), (1) 

which implies that if and only if the control qubit (the first qubit) is in the state |— ), a phase flip happens to the 
state |— ) of the target qubit (the second qubit), but nothing happens otherwise. 

The NTCP gate considered here consists of n two-qubit CP gates [Fig. 1(a)]. Each two-qubit CP gate involved in 
this NTCP gate has a shared control qubit (labelled by 1) but a different target qubit (labelled by 2, 3, or n + 1). 
According to the transformation (1) for a two-qubit CP gate, it can be seen that this NTCP gate with one control 
qubit (qubit 1) and n target qubits (qubits 2,3,...,n-|-l) can be described by the following unitary operator 



n+l 

C/=na.-2|-i-,)(-i-,|), (2) 

where the subscript 1 represents the control qubit 1, while j represents the jth target qubit, and Ij is the identity 

operator for the qubit pair (1, j), which is given by Ij = ^^.^ |riSj) (riSj|, with r, s €E {+, — }. One can sec that the 
operator (2) induces a phase flip (from the -|- sign to the — sign) to the logical state |— ) of each target qubit when 
the control qubit 1 is initially in the state |— ), and nothing happens otherwise. 

It should be mentioned that the NTCP gate can be defined using the eigenstates |0) and |1) of the Pauli operator 
(Tz- However, we note that to construct a n-target controlled-NOT gate (based on the NTCP gate defined in the 
eigenstates of dz), a Hadamard gate acting on each target qubit before and after the NTCP gate (i.e., a total of 2n 
Hadamard gates) would be required. In contrast, the construction of the n-target controlled-NOT gate, using the 
NTCP gate defined in the eigenstates of Ux, requires only 2 Hadamard gates [see Fig. 1(b)]. Therefore, the NTCP 
here, defined in the eigenstates of ax, makes the procedure for constructing a n-target controlled-NOT gate much 
simpler. 



III. MODEL AND UNITARY EVOLUTION 



Consider (n + l) qubits interacting with the cavity mode and driven by a classical pulse. In the rotating- wave 
approximation, the Hamiltonian for the whole system is 



with 



H = Hq + Hi + H2 , 



(3) 



Ho 

Hi 
H2 



-Sz + HuJcCi'a, 



hg {aS+ + a''S-). 



(4) 

(5) 
(6) 
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The Hamiltonian (3), together with the Hamihonians (4-6), can be implemented when the qubits are atoms [40,41], 
quantum dots or superconducting devices (e.g., see section V below). Here, Hq is the free Hamiltonian of the qubits 
and the cavity mode, Hi is the interaction Hamiltonian between the qubits and the classical pulse, and H2 is the 
interaction Hamiltonian between the qubits and the cavity mode. In addition, wq is the transition frequency between 
the two levels |0) and |1) of each qubit; a (a^) is the photon annihilation (creation) operator of the cavity mode with 
frequency lJc] 9 is the coupling constant between the cavity mode and each qubit; f2, w, and tp are the Rabi frequency, 
the frequency, and the initial phase of the pulse, respectively; and Sz, S-, and S+ are the collective operators for the 
(n + 1) qubits, which are given by 

n+l n+1 n+1 

j=i i=i i=i 

where az,j = \0j) {Oj\ — , cr~ = |0j) {lj\ , and cr^ = {Oj\ , with \0j) and {j = 1, 2, being the 

ground state and excited level of the jth qubit. In the interaction picture with respect to Hq, the above Hamiltonians 
Hi and H2 are rewritten respectively as (assuming co = lvq) 

Hi = '^(e''Ps_+e-'^S+), (8) 

H2 = hg{e'^*aS+ + e-'^*a'<S-) , (9) 

where 

S = u>o-uJc (10) 

is the detuning between the |0) -H- |1) transition frequency luq of each qubit and the frequency lJc of the cavity mode. 

We now consider two special cases: ip = n and negative detuning (5 < 0, as well as tp = and positive detuning 
5 > 0. The detailed discussion of these two cases will be given in the following subsections HI. A and III.B. The results 
from the unitary evolution, obtained for these two special cases, will be employed by the two alternative approaches 
(presented in section IV) for the gate implementation. 

A. Case for pulse phase Lp = tt and detuning S < 

In this subsection, we consider the negative detuning case 6 < [Fig. 2(a)]. When (f = n, the above Hamiltonian 
(8) reduces to 

Hi = -^S,, (11) 

where 

n+1 

Sx = S- + S+ = ^ ax J, (12) 
with ax,j = (^J + cr/- Performing the unitary transformation \tjj (t)) = e~*^^* Itjj' (t)) , we obtain 

,nd_Km^H,\^'^t)), (13) 



with 



H2 = ex.p[iHit/n]H2ex.p[-iHit/h] 
2 



i {Sz -S.+ S+) e-^"* 



-i(5, + 5_-5+)e^"* 



+ H.C., (14) 



Assuming that SI ^ \6\ ,g, we can neglect the fast-oscillating terms [40-42]. Then the Hamiltonian (14) reduces to 
[40-42] 

:ff2 = §(e'^*a + e-^^V)5.. (15) 
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FIG. 2: Illustration of different detunings between the cavity mode frequency Uc and the (lubit transition frequency ljq. Here, 
Wo, Wc, and w are the qubit transition frequency, the cavity mode frequency, and the pulse frequency, respectively. In (a), the 
detuning 5 = too — ujc < 0. In (b), the detuning S' = ioo — coc > 0. To implement an NTCP gate, we propose two alternative 

approaches (sec section IV), each one with three steps. The first stop for cither method requires a negative detuning {S < 0), 
as shown in (a). For either method, the second step requires a positive detuning {5' > 0), as shown in (b). 



The evolution operator for the Hamiltonian (15) can be written in the form [19,40,43,44] 

u (t) = e-^^W^^e-*^«^^'^e-*^*(*)^^'^', (16) 

where 

E(0^|£e-V = J^(e--l), 



A{t) = i f B{t') '^^]P dt' 
Jo 



9l 
4^ 



t+l(e---l 



dt' 

(17) 



i5 

When the evolution time t satisfies t = r = 27r/ |^| , we have B{t) =Q and A (r) = g'^r/ (45). Then we obtain 

M(r)=exp(aT52), (18) 

where A = —g^/ (45) > 0. The evolution operator of the system (in the interaction picture with respect to Hq) is thus 
given by 

U (r) = e-'^^^l^u (r) = g'^^^^/^giArSj _ ^ 

In section IV, we propose two alternative approaches for the NTCP gate implementation, each one involving three 
steps. The evolution operator U in Eq. (19) will be needed for the first step for either one of the two alternative 
approaches. 

B. Case for pulse phase <^ = and detuning 5 > 

In the following, we consider the positive detuning case 5 > 0. The detuning 5 can be adjusted from 5 < to 

5 > by changing either the qubit transition freqiiency or the cavity mode frequency w^j or both. In general, the 
qubit-cavity coupling constant varies when the detuning changes. To distinguish this case (^p — ^ and 5 > 0) from the 
previous case {(p = w and S < 0), we replace the previous notation fi, S, and g by fl' , 6', and g', respectively [see Fig. 
2(b)]. For simplicity, we use the same symbols lo, ljq, and cOc for the pulse frequency, the qubit transition frequency, 
and the cavity mode frequency [Fig. 2(b)]. 

Suppose now that qubit 1 is largely detuned (decoupled) from both the cavity mode and the pulse. This can be 
achieved by adjusting the level spacing of qubit 1 (this tunability of energy levels can be achieved in different types of 
qubits, e.g., solid-state qubits [30-34]), or by moving qubit 1 out of the cavity mode (e.g., trapped atoms [20,45,46]). 
Thus, after dropping the terms related to qubit 1 (i.e., the terms corresponding to the index j = 1) from the collective 
operators S+ = J2]=i '^t ^^'^ ^~ ~ Sj=i '^J Hamiltonians (8) and (9), and replacing fi, 5, and g (for the case of 
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ip ~ IT and 5 < 0) with f2', 5' , and g' (for the case of <^ = and 5 > 0), we obtain from the Hamiltonians (8) and (9) 
(assuming ui = ojq) 

H[ = !}^[e^'^S'_ + e-^^S'^), (20) 

H'^ = %' (e'^'*a5; + e-'*'*a+5L) , (21) 

which are written in the interaction picture with respect to Hq in Eq. (4). Here, S'_ and S'^ are the coUective operators 
for the n qubits (2, 3, n + 1), which are given by 



n+1 n+1 

^- = E^7' ^;=E4- (22) 

j=2 j=2 

In addition, the detuning 6' is 

S' =uio-LJc> 0. (23) 
When choosing <^ = 0, the Hamiltonian (20) reduces to 

H'l = ^S'^, (24) 

where 

n+1 

5; = y +5; = ^a,,,-. (25) 

Note that the Hamiltonian (24) has a form similar to Eq. (11) and that the Hamiltonian (21) has a form similar to 
Eq. (9). Therefore, it is straightforward to show that under the condition il' ^ 5', g' , when t = t' = 2Tr/6', the 
evolution operator of the qubits (in the interaction picture with respect to Hq) is 

U' (r') = eM-i^'r'S'j2) expi-iX'r' S'^^), (26) 

where A' = 5'^/ {4:6') > 0. The evolution operator U' in Eq. (26) will be needed for the second step of either one of 
the two alternative approaches (presented in section IV below) for implementing the NTCP gate. 

One can see that the operator described by Eq. (19) or Eq. (26) does not include the photon operator a or of 
the cavity mode. Therefore, the cavity mode can be initially in an arbitrary state (e.g., in a vacuum state, a Fock 
state, a coherent state, or even a thermal state). 

For the gate realization, we will also need to have the qubits decoupled from the cavity mode and apply a resonant 
pulse to each qubit. Suppose that the Rabi frequency for the pulse applied to qubit 1 is fii while the Rabi frequency 
for the pulses applied to qubits (2, 3, n + 1) is Sir- The initial phase for each pulse is <^ = 0. The free Hamiltonian 
of the qubits and the cavity mode is the Hq given in Eq. (4). In the interaction picture with respect to Hq, the 
interaction Hamiltonian for the qubit system and the pulses is given by 

For an evolution time r, the time evolution operator for the Hamiltonian (27) is 

U (r) = exp exp [-inrTS'^/2] . (28) 

The evolution operator U in Eq. (28) will be needed for the third step for either one of the two alternative 
approaches below. 

IV. IMPLEMENTATION OF THE NTCP GATE 



The goal of this section is to demonstrate how the NTCP gate can be realized based on the unitary operators (19), 
(26), and (28) obtained above. We will present two alternative methods for the gate implementation: one method 
based on the adjustment of the qubit transition frequency and another method which works mainly via the adjustment 
of the cavity mode frequency. 



7 



Qubiti 



Qubits (2,3,.. .,n+l) 



(a) 



Small detuning S 



(a') 



Small detuning S 



vwwv^ 

n 

|o)- 



<5<0 



^ = ^r 

n 



(b) 



Large detuning A 

10 



(b') Small detuning d' 



1 



(i) = 

n' 

|o)- 



(c) 



Large detuning A 

|i>' 



CO =Q)(j COq 
VWWW 

n, 



|o>- 



(C) 



Large detuning A 



(S = 

CO = (0q a)Q 

VWWV». 

|o)_ 



FIG. 3: Change of the qubit transition frequency too (or the qubit level spacings) of each qubit during a three-step NTCP gate. 
The three figures on the left side corresponds to the control qubit 1 while the three figures on the right side corresponds to 
each of the target qubits (2, 3, n + 1). Figures (a) and (a'), figures (b) and (b'), and figures (c) and (e') correspond to the 
operations of step (i), step (ii), and step (iii), respectively. In each figure, the two horizontal solid lines represent the qubit 
energy levels for the states jO) and |1) ; cjq is the qubit |0) o jl) transition frequency; oJc is the cavity mode frequency; w is the 
pulse frequency; ip is the initial phase of the pulse; and Q, il' , ili, and fir are the Rabi frequencies of various applied pulses. In 
addition, 5 and 5' are the small detunings of the cavity mode with the |0) |1) transition, which are given hy d = uo — coc < 
and 5' = Wo — Wc > 0; while A = wo — Wc represents the large detuning of the cavity mode with the |0) |1) transition. Note 
that the cavity mode frequency ujc is kept fixed during the entire process, but the qubit transition frequency ojq is adjusted to 
achieve a different detuning S, 5' , or A for each step. 



A. NTCP gate via adjusting the qubit transition frequency 

Let us consider (n + 1) qubits placed in a single-mode cavity or coupled to a resonator. For this first method, the 
cavity mode frequency cOc is kept fixed. The operations for the gate implementation, and the unitary evolutions after 
each step of operation, are listed below: 

Step (i): Apply a resonant pulse (with = tt) to each qubit. The pulse Rabi frequency is fl. The cavity mode is 
coupled to each qubit with a detuning S < [Fig. 3(a) and Fig. 3(a')]. This is the case discussed in subsection III. A. 
Thus, the U in Eq. (19) is the evolution operator for the qubit system for an interaction time r = —2tt/6. 

Step (ii): Adjust the qubit transition frequency for qubits (2, 3, n+ 1), such that the cavity mode is coupled to 
qubits (2, 3, n + 1) with a detuning 5' > [Fig. 3(b')]. Apply a resonant pulse (with ip = 0) to each of the target 
qubits (2, 3,...,n+ 1) [Fig. 3(b')]. The pulse Rabi frequency is now f2'. In addition, adjust the transition frequency 
of qubit 1 such that qubit 1 is decoupled from the cavity mode and the pulses applied to qubits (2, 3, n + 1) [Fig. 
3(b)]. It can be seen that this is the case discussed in subsection III.B. Thus, the U' in Eq. (26) is the evolution 
operator for the qubit system for an interaction time r' = 2Tr/6' . 
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The combined time evolution operator, after the above two steps, is 



U{t + t') = U'{t')U{t) 

^ g-iji'T's;/2g- 

= exp [i^T {S, - S',) /2] exp [iXt {Sl - S'^)] . 



(29) 



In the last line of Eq. (29), we assumed fir = VI't' (i.e., — 0/5 = U,' /8') and At = A't' (i.e., —g/5 = g' /S'), which 
can be achieved by adjusting the detunings 6 and S' (i.e., changing the qubit transition frequency) as well as the 
Rabi frequencies il and O' (i.e., changing the intensity /amplitude of the pulses). Note that Sx — S'^ = ax,i and 
— S'^ = I + 2ax,iS'^, where I is the identity operator for qubit 1. Thus, Eq. (29) can be written as 

/7 (r + r') = exp [inTax^i/2] exp [i2XTaxsS'J . (30) 

Step: (iii) Leave the transition frequency unchanged for qubit 1 [Fig. 3(c)] but adjust the transition frequency 
of qubits (2, 3,...,n+ 1) [Fig. 3(c')], such that the cavity mode is largely detuned (decoupled) from each qubit. In 
addition, apply a resonant pulse (with (p ~ 0) to each qubit. The Rabi frequency for the pulse applied to qubit 1 is 
now Oi [Fig. 3(c)] while the Rabi frequency of the pulse applied to qubits (2,3,...,n+ 1) is il^ [Fig. 3(c')]. In this 
case, the interaction Hamiltonian for the qubit system and the pulses is given by Eq. (27) and thus the time evolution 
operator is the U in Eq. (28) for an evolution time t given above. 

It can be seen that after the above three-step operation, the joint time-evolution operator of the qubit system is 



C/(2t + t') = U{t)U{t + t') 

g-i(ai-a)r<7x,i/2g-iQ^rSi/2gi2Ar<Ta,,iSi 

(31) 



Under the following condition 

Oi = 4An + n = -ng'^/S + O, 

= AX = -g'^/5, (32) 

(which can be obtained by adjusting the Rabi frequencies fi, Oi and fl^), Eq. (31) can be then written as 

n+l 

i7(2r + r')= (I'i)' (33) 

with 

Up (1, j) = exp [-i2AT {ax,i + <Jx,j - (Jx,io-x,j)] ■ (34) 
It can be easily shown that for the qubit pair (1, j), we have 

uAi,j)\+i)\+j) = 

Up{l,j)\-i)\+j) = 

Up{l,j)\-i)\-j) = exp(i8Ar)|-i)|-,), (35) 

where an overall phase factor cxp(— i2Ar) is omitted. Here and below, |-|-i) = (|0i) + |li)) /a/2, and | — i) = 
(|0i) — |li)) /V2, are the two eigenstates of the Pauli operator ax,i for qubit 1; while — {\0j) + /V2, 

and |— j) = {\0j) — /v^, are the two eigenstates of the Pauli operator axj for qubit j {j = 2, 3, or n + 1). By 
setting 8At = {2k + 1) tt, i.e., 

457^2 = 2k + l (36) 

{k is an integer), we obtain from Eq. (35) 

«7p(l,j)|+i)|+i) = 

c/p(i,i)l+i)|-,) = 
c^p(i,j)l-i)l+.-) = 

UAlJ)\-i}\-j} = (37) 
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which shows that a quantum phase gate described by 

Up{l,j) = Ij-2\-i-j){-i-j\, (38) 

is achieved for the qubit pair Here, Ij is the identity operator for the qubit pair which is given by 

Ij = \kilj) {kiljl with k. I e {+, — } . Note that the condition (36) can be achieved by adjusting the detuning 5 
(i.e., via changing the qubit transition frequency ojq). 
Combining Eqs. (33) and (38), we finally obtain 

n+l 

U (2r + r') = n - 2 Hi-) ' (39) 

i=2 

which demonstrates that n two-qubit CP gates are simultaneously performed on the qubit pairs (1,2), (1,3),..., and 
(l,n + 1), respectively. Note that each qubit pair contains the same control qubit (qubit 1) and a different target 
qubit (cither qubit 2,3. ....or n + l). Hence, an NTCP gate with n target qubits (2, 3, n + 1) and one control qubit 
(qubit 1) is obtained after the above three-step process. 

From the description above, one can see that the method presented here has an advantage: it does not require the 
adjustment of the cavity-mode frequency. 

B. NTCP gate via mainly adjusting the cavity-mode frequency 

In this subsection, we present a different way for realizing the NTCP gate, which is mainly based on the adjustment 
of the cavity mode frequency lOc- 

Let us consider (n + l) qubits placed in a single- mode cavity or coupled to a resonator. Note that now the transition 

frequency ojq of the n target qubits (2, 3, n+l) is kept fixed during the following entire process. The operations 
for the gate implementation and the unitary evolutions after each step are listed below: 

Step (i): Apply a resonant pulse (with ip = tt) to each qubit. The pulse Rabi frequency is ft. The cavity mode is 
coupled to each qubit with a detuning S < [Fig. 4(a) and Fig. 4(a')]. It can be seen that this is the case discussed 
in subsection HI. A. Thus, the U in Eq. (19) is the evolution operator for the qubit system for an interaction time 
r = -2Tr/S. 

Step (ii): Leave the transition frequency coq of qubits (2, 3,...,n+ 1) unchanged while adjusting the cavity mode 
frequency uic- The cavity mode is coupled to qubits (2, 3,...,n+ 1) with a detuning S' > [Fig. 4(b')]. Apply a 
resonant pulse (with = 0) to each of qubits (2,3,...,n+ 1) [Fig. 4(b')]. The pulse Rabi frequency is now ft'. In 
addition, adjust the transition frequency of qubit 1, such that qubit 1 is largely detuned (decoupled) from the cavity 
mode as well as the pulses applied to qubits (2,3, ...,n + 1) [Fig. 4(b)]. One can see that this is the case discussed 
in subsection III.B. Hence, the U' in Eq. (26) is the evolution operator for the qubit system for an interaction time 
r' = 27r/S'. 

Step (iii): Leave the transition frequency of qubits (2, 3, ...,n+ 1) unchanged while adjusting the transition 

frequency of the control qubit 1 back to its original setting in step (i) [Fig. 4(c)]. Adjust the cavity mode frequency 
uJc, such that the cavity mode is largely detuned (decoupled) from each qubit [Fig. 4(c) and Fig. 4(c')]. In addition, 
apply a resonant pulse (with (f = 0) to each qubit. The Rabi frequency for the pulse applied to qubit 1 is now fli 
[Fig. 4(c)] while the Rabi frequency of the pulses applied to qubits (2,3, ...,n + 1) is now ilr [Fig- 4(c')]. Therefore, 
the interaction Hamiltonian for the qubits and the pulses is Hi in Eq. (27) and thus the time evolution operator is 
the U in Eq. (28) for an evolution time r given above. 

Note that the time evolution operators U, U' , and U, obtained from each step discussed in this subsection, are the 
same as those obtained from each step of operations in the previous subsection. Hence, it is clear that the NTCP 
gate can be implemented after this three-step process. 

From the description above, it can be seen that the transition frequency uq for each one of the n target qubits 
(2, 3,...,n+ 1) remains unchanged during the entire operation. Thus, adjusting the level spacings for the n target 
qubits (2, 3, n+l) is not required by the method presented in this subsection. What one needs to do is to adjust 
the cavity mode frequency and the level spacings of the control qubit 1. 

C. Discussion 

We have presented two alternative methods for implementing an NTCP gate. Note that when coupled to a cavity 
(or resonator) mode and driven by a classical pulse, physical qubit systems (such as atoms, quantum dots, and 
superconducting qubits) have the same type of interaction described by the Hamiltonian (3) or a Hamiltonian having 
a similar form to Eq. (3), from which the four Hamiltonians (8), (9), (20), and (21), i.e., the key elements for the 
proposed NTCP gate implementation, are available. Therefore, the two alternative methods above are quite general. 
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FIG. 4: Change of the cavity mode frequency ujc and the transition frequency wo (or the level spacings) of qubit 1 during a 
three-step NTCP gate. This second method is an alternative way to implement the NTCP gate, which is different from the 
steps shown in Fig. 3. Note that the transition frequency cjq for qubits (2, 3, n + 1) remains unchanged during this method. 
The three figures on the left side correspond to the control qubit 1, which (from top to bottom) are for the operations of step (i), 
step (ii), and step (iii), respectively. The three figures on the right side correspond to each of the qubits (2, 3, n + 1), which 
(from top to bottom) are respectively for the operations of step (i), step (ii), and step (iii). In each figure, the two horizontal 
solid lines represent the qubit levels \0) and |1) ; i^o is the qubit transition frequency; lu^ is the cavity mode frequency; cj is the 
pulse frequency; ip is the initial phase of the pulse; and f2, Q,' , ili, and ilr are the Rabi frequencies of the pulses (for various 
steps). In addition, S and S' are the small dctunings of the cavity mode with the |0) -f-^- jl) transition, which arc given by 
6 = u)o — ujc < and 5' = u)o — u)c > 0; while A = ujq — ujc represents the large detuning of the cavity mode with the |0) <r^ |1) 
transition. 



which can be applied to implementing the NTCP gate with superconducting qubits, quantum dots, or trapped atomic 
qubits. 

As shown above, the first method is based on the adjustment of the qubit level spacings while the second method 

works mainly via the adjustment of the cavity mode frequency. For solid-state qubits, the qubit level spacings can 
be rapidly adjusted (e.g., in 1 ~ 2 nanosecond timescale for superconducting qubits [47]). And, for superconducting 
resonators, the resonator frequency can be fast tuned (e.g., in less than a few nanoseconds for a superconducting 
transmission line resonator [36]). 

In the appendix, we provide guidelines on protecting qubits from leaking out of the computational subspace in the 
presence of more qubit levels. As discussed there, as long as the large detuning of the cavity mode with the transition 
between the irrelevant energy levels can be met, the leakage out of the computational subspace induced due to the 
cavity mode can be suppressed. In addition, we should mention that since the applied pulses are resonant with the 
]0) |1) transition, the coupling of the qubit level ]0) or ]1) with other levels, induced by the pulses, is negligibly 
small. Therefore, the leak-out of the computation subspace, due to the application of the pulses, can be neglected. 

As discussed above, during the gate operation for either of the two methods, decoupling of the qubits from the 
cavity was made by adjusting the qubit frequency luq or the cavity frequency lUc- This particularly applies to solid 
state qubits because the locations of solid-state qubits (such as superconducting qubits and quantum dots) are fixed 
once they are built in a cavity or resonator. However, it is noted that for trapped atomic qubits, decoupling of qubits 
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from the cavity can be made by just moving atoms out of the cavity and thus adjustment of the level spacings of 
atoms or the cavity frequency is not needed to have atoms decoupled (largely detuned) from the cavity mode (e.g., 
as shown in Section VI). 

For both step (i) and step (ii) of the two methods above, it was assumed that the detuning 6 {6') is much smaller 
than the qubit Rabi frequency fl {fl') and that the Rabi frequencies are equal for all the qubits. In reality, the Rabi 
frequency for each qubit might not be identical, and the effect of the small Rabi frequency deviations on the gate 
performance may be magnified when the operation time is much longer than the inverse Rabi frequency. Hence, to 
improve the gate performance, the Rabi frequency deviations should be as small as possible, which could be achieved 
by adjusting the intensity of the pulses applied to the qubits. 

We should mention that when the mean photon number n of the cavity field satisfies n < 1, the multi- photon 
process can be neglected. In the case that the condition n < 1 does not meet, one can adjust the qubit level spacings 
or choose the qubit level structures appropriately such that the multi-photon process is negligible in the process of 
quantum operations. 

Before closing this section, it should be mentioned that the present method is based on an effective Hamiltonian 

ra+l ri+1 

Hes = ^-f^y = 2?iA ^ {ax,i + a^j - cr^,i<Tx,j) , (40) 

which can be found from Eqs. (33) and (34). One can see that this Hamiltonian contains the interaction terms between 
the control qubit (qubit 1) and each target qubit, but does not inchide the interaction terms between any two target 
qubits. Note that each term Hij in Eq. (40) acts on a different target qubit, with the same control qubit, and that any 
two terms Hij for different j's commute with each other. Therefore, the n two-qubit controUed-phase gates forming 
the NTCP gate can be simultaneously performed on the qubit pairs (1,2), (1,3),..., and {l,n + 1), respectively. 

V. REALIZING THE NTCP GATE WITH SUPERCONDUCTING QUBITS COUPLED TO A 

RESONATOR 

The methods presented above for implementing the NTCP gate are based on the four Hamiltonians (8), (9), (20) 
and (21). In this section, we show how these Hamiltonians can be obtained for the superconducting charge qubits 
coupled to a resonator. We will then show how to apply the first method above to implement the NTCP gate with 
(n + 1) charge qubits selected from A'' charge qubits coupled to a resonator (1 < n < N). A discussion on the 
experimental feasibility will be given later. 

A. Hamiltonians 

The superconducting charge qubit considered here, as shown in Fig. 5(a), consists of a small superconducting box 
with excess Cooper-pair charges, connected to a symmetric superconducting quantum interference device (SQUID) 
with capacitance Cjo and Josephson coupling energy Ejq. In the charge regime A ^ ^ Ejq » ksT 
(here, fcg. A, Ec, and T are the Boltzmann constant, superconducting energy gap, charging energy, and temperature, 
respectively), only two charge states, n = and n = 1, are important for the dynamics of the system, and thus this 
device [e.g., 30-32] behaves as a two-level system. Here, we denote the two charge states n = and n = 1 as the two 
eigenstates |+) and |— ) of the spin operator ct^. The reason for not using the eigenstates |0) and |1) of to represent 
the two charge states is that: in order to obtain the four Hamiltonians (8), (9), (20), and (21), in the following we 
will need to perform a basis transformation from the ct^ basis {|-|-) , |— )} to the cr^ basis {|0) , |1)} . The Hamiltonian 
describing the qubit is given by [30,31] 

Hg^-2E,{l-2ng)a,-Ej{^)a,, (41) 

where Ug = CgFg/ (2e), E, = / {2Cg + 4Cjo), ($) = 2£;jo cos (7r$/$o), = |+) (+| - h) (-| and = 
1+) (^1 + l~) (+1 • Here, Cg is the gate capacitance, Vg is the gate voltage, $ is the external magnetic flux applied to 
the SQUID loop, $o = h/2e is the flux quantum, and Ej ($) is the effective Josephson coupling energy. We assume 
that Vg = Vg'' + Vg'' + T^gi", where Vg'' {Vg") is the dc (ac) part of the gate voltage and V^i" is the quantum part of 
the gate voltage, which is caused by the electric field of the resonator mode when the qubit is coupled to a resonator. 
Correspondingly, we have 

ng=n''g'' + n'^ + n'^\ (42) 

where nf = CgVg^"/ (2e) , = CgU/7 (2e) , and nl"" = CgV^""/ (2e) . By inserting Eq. (42) into Eq. (41), we obtain 
the following Hamiltonian for the qubit-cavity system 

H = E^d^ - Ej ($) + hwca^a 

+AE,nl"a,+AE,n'^''a,, (43) 
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FIG. 5: (a) (Color online) Schematic diagram of a superconducting charge qubit. Here, Ejo is the Josephson coupling energy, 
Cjo is the Josephson capacitance, Cg is the gate capacitance, Vg is the gate voltage, and $ is the external magnetic flux applied 
to the SQUID loop through a control line (not shown), (b) A'^ superconducting qubits (red squares) are placed into a (grey) 
quasi-one dimensional transmission line resonator. Each qubit is placed at an antinode of the electric field, yielding a strong 
coupling between the qubit and the resonator mode. The two blue curves represent the standing-wave electric field, along the 
j/-direction. A subset of qubits from the A*' qubits, selected for the gate, are coupled to each other via the resonator mode; 
while the remaining qubits, which are not controlled by the gate, are decoupled from the resonator mode by setting their 
$ = $0/2, Vg'^ = e/Cg, and Vj"^ = to have their free Hamiltonian equal to zero. 



where = -2Ec{l - 2nf). When V^" = Vb cos (wt -t- (^) and T/i" = Vg"*" (a-f-at), the Hamiltonian (43) becomes 
[11,48] 

H = E^az - Ej ($) CT^ -I- hujco^a 

-{-Ml cos {uit + Lp)'(Tz + hg {a + a)^ CTz, (44) 

where = 2EcCgVQ/ {he) is the Rabi frequency of the ac gate vohage and g = 2EcCgVQ^ / (he) is the coupling 
constant between the charge qubit and the resonator mode. 

Let us now consider N identical charge qubits coupled to a single-mode resonator [Fig. 5(b)]. One can select a 
subset of qubits for the gate, while the remaining qubits, which are not controlled by the gate, are decoupled from the 
resonator mode by setting their $ = 3>o/2, V^'*'^ = e/Cg, and Vg'^ = to have their free Hamiltonian equal to zero. 
Without loss of generality, we assume that the set of qubits selected for the gate are the (n -I- 1) qubits labelled by 
1, 2, and n + 1 (here, 1 < n < N). From the discussion above, it can be seen that the Hamiltonian for the (n + 1) 
qubits and the resonator mode is 

H = eJz- Ej{<!>)S^ + huJca''a 

+hil. cos {ujt + Lp)Sz + hg {a + a'') S^, (45) 

where Sz = Yl]=i ^z,j and = Z]"=i ^xj, with a^^ = (+j| - |-j) (-j| and a^j = + (+j| . By 

setting Ez — (i.e., rig'^ = 1/2) for each qubit and defining uJa/2 = Ej ($) /h, the Hamiltonian (45) reduces to 

H = -Sr + hwrO^ a 

2 

+hflcos{ojt + (p)Sz + hg{a + a^)Sz. (46) 

Define now the new qubit basis \0j) = + /a/2 and |lj) = (|+j) — \—j}) /V2, i.e., perform a basis transfor- 
mation from the azj basis {\+j) , |— j)} to the azj basis {|0j) , for the jth qubit. Thus, in the new basis, the 
Hamiltonian (46) becomes 



H — Hq + Hi + H2 , 



(47) 
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where 

Ho = -^S, + fuv,a^a, (48) 

Hi = hSlcos{ut + (p)Sx, (49) 
H2 = hg{a + a^)S^. (50) 

Here, the collective operators ^2 and Sx are the same as those given in Eq. (7) and Eq. (12). In the interaction picture 
with respect to Ho, we obtain from Eqs. (49) and (50) (under the rotating- wave approximation and assuming co = loq) 

Hi = ^ (e*^5_ + e-^^S+) , (51) 
H2 = fig{e'^*aS+ + e-'^*a^S-) , (52) 

where the collective operators S'_ and 5+ are the same as those given in Eq. (7), and S = ojq — uic < 0- 

Note that Wo = 4£'jo cos (7r<I>/$o) /ft- Hence, the qubit transition frequency u>o can be adjusted by changing the 
external magnetic flux <I> applied to the SQUID loop of the charge qubit. 

Wc! now turn off the ac gate voltage applied to the charge qubit 1 (i.e., setting Vg'^ — for the charge qubit 1) 
and adjust the transition frequency coo of the charge qubit 1 to have qubit 1 decoupled (largely detuned) from the 
resonator mode. In this way, we can drop the terms corresponding to the index j = 1 from the collective operators 
S-^- = cTj^ and S- = X^j^i '^J involved in Hamiltonians (51) and (52). In addition, adjust either the transition 

frequency u)o of qubits (2, 3, n + 1) or the resonator frequency Wc, to achieve a detuning 6' = ujq — ujc > 0, and set 
an ac gate voltage Vg'^ = Vq cos {cut + (p) (with w = Wq) for each of qubits (2, 3, n + 1). The Rabi frequency fl' for 
each ac gate voltage (i.e., the pulse) is given by il' = ^EcCgV^j (he) . After replacing f2, 6, and g, with fi', 6' , and g', 
respectively; we can obtain from Eqs. (51) and (52) 

H[ = !^(e'^S'_+e-"^S'+), (53) 
H'2 = Kg' (e'*'*a5; + e-"'*aty ) , (54) 

which are written in the interaction picture with respect to Hq in Eq. (48). Here, the collective operators S'_ and S'_^ 
are the same as those given in Eq. (22). 

One can see that the four Hamiltonians (51-54) obtained here have the same forms as the Hamiltonians (8), (9), 
(20), and (21), respectively. Hence, the NTCP gate can be implemented with charge qubits coupled to a resonator. 
A more detailed discussion on this is given in the next subsection. 

B. NTCP gates with charge qubits coupled to a resonator 

Following the first method introduced in the previous section (see section IV. A), we now discuss how to implement 
the NTCP gate with (n + 1) charge qubits (1, 2, n + 1), coupled to a superconducting resonator. To begin with, 
it should be mentioned that: (a) for each step of the operations, the dc gate voltage V^'^'^ for each one of qubits (1, 2, 
n + 1) is set by V^'^'^ = e/Cg, such that E^^ = for each qubit; and (b) the resonator mode frequency lOc is fixed 
during the entire operation. The three-step operations for the gate realization are illustrated in Fig. 6. 

From Fig. 6, it can be seen that there is no need of adjusting the resonator mode frequency ujc. Hence, the procedure 
here for the gate realization is an extension of the first method introduced above. Note that the three time-evolution 
operators C/, C/', and U, obtained from each step, are the same as those obtained from each step in subsection IV. A. 
Therefore, following the same discussion given there, one can easily see that the NTCP gate can be implemented 
with {n + 1) charge qubits (i.e., the control charge qubit 1, as well as the n target charge qubits 2,3, and n + 1). 
Namely, after the three-step process in Fig. 6, a phase flip (i.e., |— ) — )• — |— )) on the state |— ) of each target charge 
qubit is achieved when the control charge qubit 1 is initially in the state |— ), but nothing happens to the states |-|-) 
and |— ) of each target charge qubit when the control charge qubit 1 is initially in the state |+) . 

From the above description, one can see that the spin operator cr^ is identical to Ux- In other words, the states 
|-|-) and |— ) are the eigenstates of both operators az and Gx- Thus, the NTCP gate (implemented with charge qubits 
here) is actually performed with respect to the eigenstates |-|-) and |— ) of the spin operator CTz (i.e., the two charge 
states n = and n = 1 above, for each qubit). 

According to the discussion in subsection IV.A, it can be found that to implement the NTCP gate, the following 
conditions need to be satisfled: (a) g,5 and O' ^ g' , S'; (b) fii = —ng'^/5 + fl and Sir = —g^/S] (e) —Q/S = SI' /5' 
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FIG. 6: Change of the qubit transition frequency wo and the ac gate-voltage frequency w during a three-step NTCP gate 
with charge qubits coupled to a resonator. The three figures on the left side correspond to the charge control qubit 1, which 
(from top to bottom) are for the operations of step (i), step (ii), and step (iii), respectively. The three figures on the right 
side correspond to the charge target qubits (2, 3, n + 1), (from top to bottom) for the operations of step (i), step (ii), and 
step (iii), respectively. In each figure, the two horizontal solid lines represent the qubit levels |0) and loc is the resonator 
mode frequency; wo is the qubit transition frequency; tj is the frequency of the ac gate voltage (i.e., the pulse); and Q. (Vb), 
O' (Vq) , f2i (Vo^), or fir (ViJ) is the function of the amplitude Vb, Vj, Vq^, or VJ" of the ac gate voltage, which can be adjusted 
by changing the ac gate- voltage amplitude; and each circle with a symbol ~ represents an ac gate voltage. In (b), the ac gate 
voltage for qubit 1 is set to zero (i.e., = 0). In addition, 5 and 5' are the small detunings of the resonator mode with the 
|0) O |1) transition, which are given by 5 = wo — Wc < and 5' = wo — Wc > 0; while A = wo — Wc represents the large detuning 
of the resonator mode with the |0) <-> |1) transition. Note that the resonator mode frequency Wc is kept fixed during the entire 
operation, but the qubit transition frequency wo is adjusted to achieve a different detuning 5,5', or A for each step. 



and —gjS = g' /S'; and (d) Ag^ /S'^ = 2k + 1. These conditions can in principle be realized because: (i) The Rabi 
frequencies V, (Vo) , ^' {Vq) , fii (Vq^) , and Qr {Vq) are respectively functions of the amplitudes Vq, Vq, Vq, and Vq of 
the ac gate voltages, which can be adjusted by changing the amplitudes of the ac gate voltages; and (ii) The detunings 
S and 5' can be adjusted by changing the qubit transition frequency ujq- 

Note that = 1/2 (i.e., — 0) was set for each qubit during the entire operation. We now give some discussion 
on the deviation from the degeneracy point n^*^ = 1/2 for each one of the qubits (1, 2, n + 1), during each step 
shown in Fig. (6). 

(a) For each one of the qubits (1, 2, n -|- 1) in step (i), we have, from Eq. (42), that 

ng = nf + cos (wi -|- 1^) -|- nc["(a + a^), (55) 

where 

< = CgVo/ (2e) = hn/ {4E,) , 
nr = CgV^V{'^e) = hg/{AE,). 

Therefore, the maximal deviation from the degeneracy point for each one of the qubits (1, 2, n-|-l)in step (i) is 

So = \nf + nl" + ng" - 1/2| =h{Vt + g)/ {AE^) . (56) 
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FIG. 7: (Color online) Proposed setup for three qubits (red squares) and a (grey) standing-wave quasi-one dimensional coplanar 
waveguide cavity (not drawn to scale). Each qubit is placed at an antinode of the electric field. The two blue curves represent 
the standing-wave electric field, along the y-direction. Vg is the gate voltage, D is the distance between any two nearest SQUID 
loops, and d is the linear dimension of each SQUID loop. 

(b) Similarly, one can find that the maximal deviation from the degeneracy point for qubits (2, 3, n -I- 1) in step 
(ii) is 

ei^n{n' + g')/{4E,). (57) 

Note that the deviation from the degeneracy point for qubit 1 is smaller than ei since Vg'^ = and thus Uq'^ — for 
qubit 1. 

(d) For step (iii), we have ng'^ — h£li/ (-iEc) for qubit 1, ng'^ = hftr/ (4£'c) for qubits (2, 3, n -|- 1), and n^}" — 0. 
Thus, it is easy to see that the maximal deviation from the degeneracy point for qubit 1 is 

e2^hni/iAE,), (58) 

and the deviation from the degeneracy point for qubits (2, 3, n -f 1) is 

£3 = hnr/ i^E,) . (59) 

In the above, we have discussed how to realize the NTCP gate with superconducting charge qubits coupled to a 
resonator. Note that the proposal here can implement multi-qubit gates, while previous proposals (e.g. ,[49-51]) using 
superconducting qubits are limited to two-qubit gates. 

C. Possible experimental implementation 

In this section we discuss some issues which are relevant for future experimental implementation of our proposal. 
For the method to work: (a) The conditions for the Rabi frequencies fi, fi', fii, and fir, which were discussed above, 
need to be met; (b) The total operation time 

iop = 2r + r' = 47r/|(5| +27r/(5' (60) 

should be much shorter than the energy relaxation time Ti and the dephasing time T2 of the qubit and the lifetime 
of the resonator mode k^^ ~ Q/^j^c, where Q is the (loaded) quality factor of the resonator; (c) The deviations £0, 
£1, £2, and £3 from the degeneracy point need to be small numbers to have the qubits working near the degeneracy 
point, such that the qubits are less affected by the low-frequency charge noises [52,53]; and (d) The direct coupling 
between SQUlDs needs to be negligible, since this interaction is not intended. It is noted that the direct interaction 
between SQUlDs can be made negligibly small as long a.s D ^ d (where D is the distance between the two nearest 
SQUlDs and d is the linear dimension of each SQUID). 

For the sake of definitiveness, let us consider the experimental feasibility of implementing a two-target-qubit con- 
trolled phase gate using superconducting charge qubits with parameters listed in Table 1 [8,54-56]. Note that in a 
recent experiment, coupling three superconducting qubits with a transmission line resonator has been demonstrated 
[57]. For a superconducting one dimensional standing- wave CPW (coplanar waveguide) transmission line resonator 
and each qubit placed at an antinode of the resonator mode (Fig. 7) , the amplitude of the quantum part of the gate 
voltage is given by [48] 

= {hj^f' [LcoT"^ , (61) 
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Charge qubit 


Gate capacitance 


Q~ 1 aF [54] 


Josephson capacitance 


CjQ ~ 300 aF [541 


Charging energy 


EJh~-il GHz [541 


Tn';prtfi*;nn rniinlin(y pTiprtrv 


t,jQi n ~j vj-Hz ^-'^J 


Energy relaxation time 


Ti ~ 7.3//S [55] 


Dephasing time 


T2 ~ 500 ns [55] 


Rabi frequency 


np/2;r~300MHz[56] 


Tuning range of the qubit 
frequency 


Avo ~350MHzf*^ 



TABLE 1: Possible experimental parameters of a charge qubit [8,54-56]. Here, flp is the typical Rabi frequency available in 
experiments, and Ai>o is the experimentally tuning range of the qubit frequency. 



Resonator 


Resonator frequency 


m,./2^ = 10 GHz 


Wave length 


Z ~10 nun 


Lengtli of tlie resonator 


L = Z 


Capacitance per unit lengtti of resonator 


c,) ~0.22aF 


Etfcctivc iclati\c dielectric constant 


f - 9 


Quality factor 


Q ~ Uf 


Tuning range of tlie resonator frequency 


Av,. ~ 740MHz 



TABLE 11: Possible experimental parameters of a resonator [36,57-59]. Here, Avc is the experimentally tuning range of the 
resonator frequency. 

where L is the length of the resonator and cq is the capacitance per unit length of the resonator. Therefore, the 
coupling constant g is given by 

g = 2E^Cg {he)-^ {hw^f^ {Lco)'^'^ , (62) 

showing that g does not depend on the detuning 5. Therefore, we have g = g', for which the above condition 
—0/(5 = 1 5' and —g/5 — g' 15' simply turns to the = $7' and —5 = 5' . For superconducting charge qubits with 
parameters given in Table 1, and a resonator with the parameters listed in Table II [36,57-59], a simple calculation 
gives g/2iT ~ 22 MHz, which is available in experiments (see, e.g., [58]). With a choice of —5 = 5' = 2g [corresponding 
to A; = in Eq. (36)], the total operation time top would be ~ 68 ns, which is much shorter than the dephasing time 
T2 and n^^ ^ 794 ns for a resonator with Q = 10"''. Note that a superconducting CPW resonator with a quality factor 
of Q > 10^ has been experimentally demonstrated [59]. 

Based on f2 = fi', fii = 4An + Vt, and fl^ = 4A (n = 2 for a two-target-qubit gate), we have (27r) ~ 330 MHz, 
fli/(27r) - 352 MHz, and 17^/ (27r) ~ 11 MHz for a choice of 17 - 15g, i.e., 17/ (27r) - 330 MHz. For the Rabi 
frequencies given here, we obtain eq = £1 ~ 2.75 x 10~^, £2 ~ 2.75 x 10"'^, and £3 ~ 8.59 x 10~^ for a qubit-cavity 
system with the above parameters. Therefore, the conditions for the qubits to work near the degeneracy point are 
well satisfied. 

Finally, for a resonator with ujc/ (27r) = 10 GHz, the wavelength of the resonator mode is A ^ 10 mm. For the 
charge qubits placed in a resonator as shown in Fig. 7, the distance between any two nearest SQUIDs is I? ^ A/2 ~ 5 
mm. Hence, the ratio D/d would be ~ 250 for d = 20/xm. Note that the dipole field generated by the current in each 
SQUID ring at a distance r^d decreases as r~^. Thus, the condition of negligible direct coupling between SQUIDs 
is well satisfied. 

Note that because of wo = 4_Ejo cos (7r$/$o) /ft, it can be found that for charge qubits with the parameters given 
in Table I, the transition frequency uq = a;o/27r of each qubit varies from z/q = GHz for $/$o = 1/2 to vq = 20 
GHz for $ = 0. Therefore, the choice of the resonator frequency above is reasonable. For the choice oi —5 = 5' = 2g 
above, the qubit transition frequency vq would be ~ 9.956 GHz for the detuning 5 = —2g while ~ 10.044 GHz for the 
detuning 5' = 2g. 

The above analysis shows that the realization of a two-target-qubit controlled phase gate is possible using super- 
conducting charge qubits and a resonator. We remark that a quantum-controlled phase gate with a larger number of 
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(c) 

FIG. 8: (Color online) Proposed set-up for an n-target-qubit control phase (NTCP) gate with (n + 1) identical neutral atoms 
and a cavity. For simplicity, only five atoms are drawn here. Each atom can be either loaded into the cavity or moved out of 
the cavity by one-dimensional translating optical lattices [20,45,46]. The atom in the middle represents atom 1 (the control 
qubit), while the remaining atoms play the role of target qubits. 

target qubits can in principle be obtained by increasing the length of the resonator since the total operation time top 
is independent of the number of target qubits n. 

VI. NTCP GATE WITH ATOMS USING ONE CAVITY 

Consider (n -I- 1) identical two- level atoms (l,2,...,n-|-l). The two levels of each atom are labelled by |0) and 
|1) . The transition frequency of each atom is denoted as ujo- Each atom is trapped in the periodic potential of a 
one-dimensional optical lattice and can be loaded into or moved out of the cavity by translating the optical lattice 
[20,45,46]. The NTCP gate can be realized using a procedure illustrated in Fig. 8. The operations shown in Fig. 8 
are as follows: 

• Figure 8(a): Move atoms (l,2,...,n-|-l) into the cavity and then apply a classical pulse (with an initial phase 
if — TT and frequency uj — ujq) to the atoms. The cavity mode is coupled to the ]0) — > ]1) transition of each atom, 
with a detuning S ^ loq — lOc < 0, which can be achieved by prior adjustment of the cavity mode frequency [35]. The 
time-evolution operator for this operation is the U in Eq. (19) for an interaction time t = —27t/S. 

• Figure 8(b): Move atoms (1,2, ...,n + 1) out of the cavity and then adjust the cavity mode frequency. After 
adjusting the cavity mode frequency, move atoms (2, 3, ...,n + 1) back into the cavity and then apply a classical pulse 
(with if = and uj = ujq) to them. The cavity mode frequency is adjusted such that the cavity mode is coupled to 
the ]0) —7- ]1) transition of atoms (2, 3, n + 1), with a detuning S' = ujo — uJc > 0. The time-evolution operator for 
this operation is the U' in Eq. (26) for an interaction time r' = 27r/(5'. 

• Figure 8(c): Move atoms (2, 3, n+1) out of the cavity. Apply a classical pulse (with (p — and uj — ojq) to the 
control atom 1 and a classical pulse (with the same initial phase and frequency) to the target atoms (2, 3,..., n + 1). 
The Rabi frequency for the pulse applied to atom 1 is f2i, while the Rabi frequency for the pulses applied to atoms 
(2, 3, ...,n + 1) is Qr- This operation is described by the operator U in Eq. (28) for an interaction time t. 

The three time-evolution operators U, U' , and U here are the same as those obtained from each step in subsection 
IV. A. Therefore, as long as the conditions given in subsection IV. A are met, an NTCP gate described by Eq. (39) 
is implemented with the (n + 1) atoms. Namely, after the operations shown in Fig. 8, n two-qubit controlled phase 
gates, each described by the unitary transformation in Eq. (1), are performed on the atom pairs (1, 2), (1, 3),..., and 
(1,71 -|- 1), respectively. Note that each pair contains the same control qubit (atom 1) but a different target qubit 
(atom 2, 3, or n -I- 1). 

As shown above, the present scheme for implementing the NTCP gate with atoms has the following features: 

(a) No adjustment of the level spacings for each atom is required during the entire operation; 

(b) Only one cavity is required; 

(c) The cavity mode can be initially in an arbitrary state; and 

(d) The operation time does not depend on the number of atoms. 
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For our scheme to work, the total operation time top = t + t' + Ta + 4rm should be much smaller than the cavity 
decay time k^^, so that the cavity dissipation is negligible. Here, Ta is the typical time required for adjusting the 
cavity mode frequency during step (ii) above, and Tm is the typical time required for moving atoms into or out of 
the cavity. In addition, the top needs to be much smaller than the energy relaxation time of level |1), such that the 
decoherence induced due to the spontaneous decay of the level |1) is negligible. In principle, these conditions can be 
satisfied by choosing a cavity with a high quality factor Q and atoms with a sufficiently long energy relaxation time. 

To investigate the experimental feasibility of this proposal, let us consider Rydberg atoms with principal quantum 
numbers 50 and 51 (respectively corresponding to the levels |0) and |1)). The |0) O |1) transition frequency is 
~ 51.1 GHz, the energy relaxation time of the level |1) is T,. ^ 3 x 10^^ s, and the coupling constant is 
5 = 27r X 50 KHz [35,60]. Now set —5 = 5' = 2g, and assume g ~ g'. With the choice of these parameters, the 
time needed for the entire operation is top ~ 65 /xs for Ta = Tm ~ 1 /xs, which is much shorter than T^.. The cavity 
mode frequency is lUc/^tt ^ 51.2 GHz for the negative qubit detuning 5 = —2g while ^ 51 GHz for the positive qubit 
detuning 5' = 2g. To estimate the lifetime of the cavity photon, we here consider the conservative case of a larger 
cavity frequency Wc/27r ~ 51.2 GHz, for which the lifetime of the cavity photon is Tc = Q/uJc ~ 622 /is for a cavity 
with Q = 2 X 10^, which is much larger than the total operation time top- Note that cavities with a high Q ~ 10^° 
have been demonstrated in experiments [61]. Thus, the present proposal might be realizable using current cavity 
QED setups. 

VII. CONCLUSIONS 

We have presented two different methods for the proposed NTCP gate implementation. The two methods are 

quite general, which can be applied to physical systems such as trapped atoms, quantum dots, and superconducting 
qubits. For the two methods, we have provided guidelines on how to protect multi-level qubits from leaking out of 
the computational subspace. 

Using a concrete example, we have shown how to apply the first method to implement the NTCP gate with super- 
conducting qubits coupled to a resonator. In addition, we have discussed the experimental feasibility of performing 
a two-target-qubit controlled phase gate with superconducting qubits coupled to a one-dimensional transmission line 
resonator. Our analysis shows that the realization of this gate is possible within current technologies. How well this 
gate would work in light of experimental errors should be further investigated elsewhere for each particular set-up or 
implementation. This is beyond the scope of this theoretical work. 

We have shown how to extend the second method to implement the NTCP gate with trapped atomic qubits by 
using one cavity. Interestingly, as shown above, there is no need to adjust the level spacings of atoms during the entire 
operation, and decoupling of atoms with the cavity can be easily achieved by just moving atoms out of the cavity. 

In summary, we have presented a general proposal to implement a NTCP gate with qubits in a cavity or coupled to 
a resonator. As shown above, the present proposal has the following features: (i) The n two-qubit CP gates involved 
in the NTCP gate can be simultaneously performed; (ii) The operation time required for the gate implementation is 
independent of the number of the target qubits, thus it does not increase with the number of qubits; (iii) The gate 
is insensitive to the initial state of the cavity mode, therefore no preparation for a specified initial state of the cavity 
mode is needed; (iv) No measurement on the qubits or the cavity mode is needed and thus the operation is simplified, 
(v) The gate realization requires only three steps of operations. 
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APPENDIX: HOW TO PROTECT QUBITS FROM LEAKING TO HIGHER ENERGY LEVELS 

Let us now discuss how to protect qubits from leaking out of the computational subspace in the presence of more 
qubit levels. Generally speaking, we need to consider two situations, which are here denoted as cases L and S. For 
case L, the level spacing between the level |1) and the level |2) (the first level above |1)) is larger than the level 
spacing between the levels jO) and |1) . Namely, for case L, {E2 — Ei) > {Ei — Eq) , where Eq, Ei, and E2 are the 
energy eigenvalues of the levels |0) , |1) , and |2) . For case S, the level spacing between the levels |1) and |2) is smaller 
than that between the levels |0) and |1) . Namely, (£"2 — Ei) < (Ei — Eq) . For solid-state qubits, case L exists for 
superconducting charge qubits and flux qubits [31]; while case S can be applied to superconducting phase qubits [33]. 
For atomic qubits, the level structures for both cases L and S are available. 

Here and below, we define A2 = {E2 — Ei) jh — lJc, A3 = (E^ — Ei) jh — lJc, and Aa = [Ea — Ei) /h — lJc, as 
the large detuning of the cavity mode with the |1) O |2) transition, the large detuning of the cavity mode with the 
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FIG. 9: Large detuning of the cavity mode with the transition between energy levels. Figure (a) corresponds to the case L, 
while figure (b) corresponds to case S. The dots above level |2) in (a) and level |3) in (b) represent other energy levels. 



|1) -O- |3) transition, and the large detuning of the cavity mode with the |1) O \a) transition, respectively. Ea is the 
energy eigenvalue of level \a). 

A. Case L: {E2 - Ei) > {Ei - Eq) 

For case L, let us consider Fig. 9(a), where the dashed line falls within the range between the levels |1) and |2) 
when the cavity mode is slightly detuned from the |0) o |1) transition. Therefore, a large detuning A2 is needed 
to avoid the transition from the level |1) to the level |2). Note also that as long as this large detuning is met, no 
transition from level |1) to the higher energy level above the level |2) happens. To explain this, let us consider an 
arbitrary level \a) above level |2) [Fig. 9(a)]. Since the detuning is larger than the detuning A2, then the large 
detuning regime for Aq is automatically satisfied when A2 is large. As a result, the transition from the level |1) to 
the level |o) will not be induced by the cavity mode. In addition, the excitation of the level |a) , induced due to its 
coupling to the level |2) via the cavity mode, docs not occur when the level |2) is unpopulated. Hence, as long as the 
large detuning A2 is satisfied, no level above |1) will be populated and therefore the leak out of the computational 
subspace is suppressed. 

B. Case S: {E2 - Ei) < {Ei - Eq) 

For case S, let us consider Fig. 9(b), where the dashed line falls within the range between the levels |2) and |3), 
when the cavity mode is slightly detuned from the |0) O |1) transition. This can be achieved, e.g., for superconducting 
qubits by appropriately choosing the device parameters and/or adjusting the external parameters to control the level 
structures. From Fig. 9(b), it can be seen that a large detuning A2 and a large detuning A3 are needed to avoid the 
transition from the level |1) to the level |2) or the level |3). Note also that as long as these two large detunings are 
met, no transition from the level |1) to the higher energy level above the level |3) occurs. To understand this, let 
us consider an arbitrary level \a) above |3) [Fig. 9(b)]. Note that the detuning A^ is larger than the detuning A3. 
Therefore!, the large detuning regime for A^ is automatically satisfied when the cavity mode is largely detuned from 
the )1) O )3) transition. As a result, the transition from level )1) to level \a) will not be induced by the cavity mode. 
In addition, the excitation of the level |a), caused due to its coupling to the level )2) or )3) via the cavity mode, will 
not happen when levels |2) and |3) are unpopulated. Hence, as long as the large detunings A2 and A3 are met, the 
levels above |1) will not be excited by the cavity mode. 

We now give a quantitative analysis on the effect of the finite detuning of the qubit frequencies with the cavity 
mode. For simplicity, we consider the case where the cavity mode is in a single photon state and the qubit is in the 
state |1) . It is estimated that the occupation probability p2 for level |2) and the occupation probability ps for level 
|3) (induced by the photon) would be on the order of 45^2/ (4512 + '^2) ^^'^ ^Sis/ (^^is + ^3) > respectively. Here, 
gi2 is the coupling constant between the cavity mode and the |1) 4-> |2) transition, while 913 is the coupling constant 
between the cavity mode and the )1) O )3) transition. With a choice of A2 = 10(?i2 and A3 = Wgis, we have p2, 
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Ps ~ 0.04, which can be further reduced by increasing the detuning A2 or A3. Therefore, the population probabihty 
of the level |2) or |3) of the qubit can be made negligible by choosing the detuning appropriately. 
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